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ABSTRACT 



We present high-resolution spectroscopic mid-infrared observations of the cir- 
cumstellar disk around the Herbig Ae star HD 97048 obtained with the VLT 
Imager and Spectrometer for the mid-InfraRed (VISIR). We conducted observa- 
tions of mid-infrared pure rotational lines of molecular hydrogen (H2) as a tracer 
of warm gas in the disk surface layers. In a previous paper, we reported the 
detection of the S(l) pure rotational line of H2 at 17.035 /xm and argued it is 
arising from the inner regions of the disk around the star. We used VISIR on the 
VLT for a more comprehensive study based on complementary observations of 
the other mid-infrared molecular transitions, namely S(2) and S(4) at 12.278 //m 
and 8.025 //m respectively, to investigate the physical properties of the molecular 
gas in the circumstellar disk around HD 97048. We do not detect neither the 
S(2) line nor the S(4) H2 line from the disk of HD 97048, but we derive upper 
limits on the integrated line fluxes which allows us to estimate an upper limit 
on the gas excitation temperature, T^x < 570 K. This limit on the temperature 
is consistent with the assumptions previously used in the analysis of the S(l) 
line, and allows us to set stronger contraints on the mass of warm gas in the 
inner regions of the disk. Indeed, we estimate the mass of warm gas to be lower 
than 0.1 Mjup. We also discuss the probable physical mechanisms which could 
be responsible of the excitation of H2 in the disk of HD 97048. 



Subject headings: stars: 
circumstellar matter — 
stars 



pre-main sequence — stars: individual (HD97048) — (stars:) 
(stars:) planetary systems: protoplanetary disks — infrared: 
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Introduction 



Disks around young stars are a natural outcome of tlie star formation process and the 
place for planet formation. At the present time, a significant effort has been done on the 
study of the dust in disks. However, the dust only represents a tiny fraction of the disk 
mass (~1%), and it is thus mandatory to deeply study the gas phase in disks in order 
to set stronger constraints on the giant planets formation process. Molecular hydrogen 
(H2) is the most abundant molecule in the circumstellar (CS) environments of young stars 
and is supposed to be the key element of giant planet formation, thus its diagnostics are 
promising. Indeed, the detection of H2 provides the most direct information about the 
gaseous content of disks, setting limits on the timescales for the dissipation of CS matter 
and possibly planet bui l ding. H? has been observ ed in CS environments at ultraviolet (e.g . 



Johns-Krull et al. 



2000 



Martin-Zaidi et al. 



20081 ) and near-infrared (e.g. 



Bary et al. 



20031 ) 



wavelengths. These observations trace hot circumstellar gas, or gas excited by fluorescent 
processes, or require specific spatial distributions for the gas to be detectable. They are 
therefore difficult to translate into gas masses. The pure rotational mid-infrared H2 lines 
are useful probes because the level populations are expected to be in local thermodynamic 
equilibrium (LTE) at the local gas temperature, and so line ratios allow determination of 
the excitation temperature and mass of the warm gas. 

HD 97048 is a nearby, relatively isolated Herbig A 0/B9 star, surrounded by a CS disk. 



19981 ). Its age 



located in the Chameleon cloud at a distance of 180 pc (Ivan den Ancker et al.l 
has been estimated from evolutionary tracks to be of the order of 3 Myrs (kindly computed 
for us by L. Testi and A. Palacio s). The VISIR {VLT Imager and Spectrometer for the 



mid-InfraRed; iLagage et al. 



2OO4J ) imaging observations of this star have revealed emission 



of PAHs (Polycyclic Ar omatic Hydrocarbo ns) at the surface of a flared disk extending 



at least up to 370 AU (ILagage et al. 



20061 ). The flaring index has been measured to be 
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1.26±0.05, which is in go od agreement with hydrostatic flared disk models (iLagage et al. 



2006 



Doucet et al. 



20071 ). This is the only Herbig star for which the flaring of the disk has 
been observed by direct imaging. This geometry implies that a large amount of gas should 
be present to support the hydrostatic structure and that the disk is at an early stage of 
evolution. This star is thus one of the best candidates to study the gas component in the 
disks of HAes. 



In a previous paper (jMartin-Zai'di et al.l 120071 . hereafter Paper I), we presented VISIR 
high-resolution spectroscopic mid-infrared observations of the circumstellar disk around the 
Herbig Ae (HAe) star HD 97048 near 17.03 fim. Although 



Carmona et al. 



(120081 ) suggested, 



from their disk model with the assumptions of LTE conditions and gas-to dust ratio of 
about 100, that mid-IR H2 lines should not be detected with the existing instruments, we 
detected the S(l) pure rotational li ne of molecular hyd rogen at 17.035 /im arising from the 



disk around HD 97048. In addition, 



Bitner et al. 



(I2OO7I ) have detected mid-IR H2 rotational 



lines in the disk of another Herbig star, namely AB Aur, using the high spectral and spatial 
resolution TEXES spectrometer. These detections demonstrate that H2 can be observed in 
the mid-IR domain when particular physical conditions exist in disks. 

Despite the fact that the line is neither spatially nor spectrally resolved, the detection 
of the S(l) line in the disk of HD 97048 revealed the presence of significant amounts of 
warm gas in the inner 35 AU of the disk. Circumstellar gas had been previously detected in 



the disk of HD 97048 by 



Acke fc van den Anckerl ( l2006l ). who showed that the [O I] emission 



arises at radii between 0.5 and 60 AU from the central star. Very recently, modeling of 
near-infrared CO emissio n showed that the inner radius of the CO emitting region is located 



at 12 AU from the star (Ivan der Plas et al. 



20081 ). These observations reinforced the claim 



that the disk of HD 97048 contains large amounts of gas both at small and large radii, as 
suggested by the flaring geometry, and that HD 97048 is a young object surrounded by a 
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disk at an early stage of evolution. Indeed photoevaporation of the gas and planet formation 
clear up the iniier pa r t of the disk in timescales expected to be lower than 3 million years 



( iTakeuchi et al. 



2005 : 



Gorti Sz HoUenbach 



2008a 



b|). The H2 S(l) line detection also implies 



that particular physical conditions for H2 are present in the inner disk surface layer. The 
analysis of these data usually assumes that the H2 excitation is in local thermodynamic 
equilibrium (LTE), and can itself thus be characterized by a single excitation temperature, 
which should be close to the g as temperature because o f the low critical densities. As shown 



by the AB Aur observations (IBitner et al. 



2007 



20081 ) ■ the H2 gas temperature could be 



significantly higher than the dust temperature in the disk surface layers of HD 97048, due 
to the dust settling or coagulation, for example. However, several competing mechanisms 
could contribute to the excitation o f molecular hydrogen such as UV and X-r ays heating. 



shocks, etc., (see review papers by 



Habart et al. 



2004 



Snow fc McCall 



20061 ) and could 



be responsible for the observed emission. We refer the reader to Paper I for a complete 
discussion about the different possible explanations. 

The detection of the other pure rotational lines of H2 is a potentially powerful tool 
to determine the excitation temperature (and thus the mass) of the warm gas, and better 
constrain the excitation of H2. In this paper, we thus present VISIR high-resolution 
spectroscopic observations of the two other pure rotational lines of H2 observable from the 
ground: the S(2) line at 12.2786 fim, and the S(4) line at 8.0250 fim. We also recall the 
main results obtained from the S(l) line detection (detailed in Paper I). 



2. Observations and data reduction 

HD 97048 was observed at 3 different epochs. The observations at 17.035 /im presented 
in the Paper I were performed in 2006 June 22, the 8.025 fim observations in 2007 April 
07, and the 12.278 fim observations in 2007 July 03. The three sets of observations were 
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obtained with the high-resolution spectroscopic mode of VISIR. The exposure time, sht 
width, and atmospheric conditions during the observations for the target star and standard 
references are presented in Tabfe [H 

For all the three observations, the standard chopping and nodding technique was used 
to suppress the large sky and telescope background dominating at mid-IR wavelengths. 
Asteroids and standard stars were observed just before and after observing HD 97048, at 
nearly the same airmass and seeing conditions as the object. In order to correct the spectra 
from the Earth's atmospheric absorption, we divided each spectrum of HD 97048 by that of 
the corresponding asteroid (which has a much better signal-to-noise ratio th an that of th e 



standard star), and used the standard stars observed and modeled spectra (ICohen et al. 



1999 1) to obtain the absolute flux calibration. In order to verify our flux calibration, 
we checked the Spitzer Space Telescope archives. HD 97048 was observed with the IRS 
spectrograph installed onboard the Spitzer Space Telescope as part of the "Cores to Discs" 
(c2d) legacy program (AOR:00056 44800, PI: Evans). The data reduction was performed 



using the c2d legacy team pipeline (ILahuis et al. 



20061) with the S15.3.0 pre-reduced (BCD) 



data. The absolute flux calibration of our VISIR data in the 3 wavelength ranges has been 
performes in order to be consistent with the Spitzer flux measurements. We estimate the 
error on the absolute flux calibration to be lower than 10% , a value which does not affect 
our analysis. 

The wavelength calibration is done by fitting the observed sky background features 
with a model of Paranal's atmospheric emission (for more details on the observation 
and data reduction techniques see Paper I). We note that A^= 0.24 mag for HD 97048 



(IValenti et al. 



20001 ) ■ thus we have not corrected the spe ctra for dust extin ction, since it is 



negligible in our wavelength range for any Av< 40 mag fiFluks et al. 



19941 ). 
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3. Results 

3.1. Non-detection of the S(2) and S(4) lines 

HD 97048 spectra show no evidence for H2 emission neither at 12.278 /im nor at 
8.025 fim (Fig. [1], middle and bottom panels). For each flux-calibrated spectrum, we 
calculated the standard deviation (cr) for wavelength ranges relatively unaffected by telluric 
absorption, and close to the wavelength of interest. The 3a upper limits on the integrated 
line fluxes were calculated by integrating over a Gaussian of FWHM equal to a spectral 
resolution element (21 kms~^) and an amplitude of about 3a flux (Fig. [T]). We assumed 
the same radial velocity for the S(2) and S(4) lines than that observed for the S(l) line and 
corrected appropriately for each epoch of observation. We thus centered the Gaussian on 
the expected wavelengths for the S(2) and S(4) lines respectively (see Table [2]). From the 
limits on integrated intensities, we estimated the upper limits on the column densities of the 
corresponding upper rotational levels of each H2 transition (see Table [2]). For this purpose. 



we flrst assumed that the line is optically thin and that the radiatio n is isotropic. 



context, the column densities are derived from the following formula fivan Dishoeck 



n this 



1992): 



lui = -r^Nu{H2)Aui ergss ^ cm sr \ (1) 
47rA 

where lui is the integrated intensity of the line, A is the wavelength of the transition 
J = u — I, Aui is the spontaneous transition probability, Nu{H2) is the column density of 
the upper rotational level of the transition. 
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3.2. Detection of the S(l) line at 17.035 fim 

We recall here the results of our analysis of the S(l) line. For more details, we refer 
the reader to the Paper I. As shown in Fig. [1] (top panel), we have detected the H2 pure 
rotational S(l) line near 17.03 fim. In the flux-calibrated spectrum, the standard deviation 
(cr) of the continuum flux was calculated in regions less influenced by telluric absorption, 
and close to the feature of interest. We deduced a 6a detection in amplitude for the line, 
corresponding to a signal-to-noise ratio of about 11a for the line, when integrating the 
signal over a resolution element (6 pixels). The line is not spectrally resolved as we can 
flt it with a Gaussian with a full width at half maximum (FWHM) equal to a spectral 
resolution element of 21 kms~^ (Fig- [I])- From our flt, we derived the integrated flux in 
the line (see Table |2]). Once the spectrum is corrected from the E arth's rotation, a nd 



Acke et al. 



2OO5I), we 



knowing the heliocentric radial velocity of HD 97048 (+21 kms~^; 
estimated, from the wavelength position of the Gaussian peak, the radial velocity of H2 to 
be about 4±2 kms^^ in the star's rest frame. We thus considered that the radial velocity 
of the H2 is compatible with zero (at the VISIR resolution) and therefore similar to that 
of the star, implying that the emitting gas is gravitationally bound to the star, and likely 
arising from the disk and not from an outflow. H2 outflows from T Tauri stars usually 
have (blueshifted) velocities of a few tens up to a hundred of kms^^ when a clear outflow 
is identifled. Otherwise, when the line is detected at zero systemic velocity and is narrow 



or unresolved as is the c ase here, it is generally attributed to the disk (e.g. 



2006 



Herczeg et al. 



Bitner et al.l 12008) . The H2 line is not resolved spatially either. Given the VISIR 



spatial resolution of about 0.427" at 17.03 fim, and the star distance (180 pc from the 
Sun), we can assess that the emitting H2 is located within the inner 35 AU of the disk. 
Assuming that the H2 gas follows the same (Keplerian) kinematics as the disk, the emitting 
gas observed with VISIR is likely not concentrated signiflcantly in the innermost AU of 
the disk (< 5 AU) otherwise rotational broadening would be observed. Indeed, near the 
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central star, the rotational velocity of the disk is of the order of a hundred of kms~^, to 
be compared to the spectral resolution of VISIR, 21 kms~^. The emitting H2 is thus more 
likely distributed in an extended region within the inner disk, between 5 AU and 35 AU of 
the disk. Assuming the emission arises from an isothermal mass of optically thin H2, we 
estimated the corresponding column densities and masses of H2 as a function of prescribed 
temperatures (for 150 K, 300 K and 1000 K; see Paper I for details about the method). In 
Table [2] are reported the integrated flux of the S(l) line, its intensity, and the corresponding 
column density for the J = 3 rotational level. Since we have a signal-to-noise ratio of about 
llcr for the line, we can deduce 9% error bars on theses values (see Table [2]). 



The estimates of the column densities of the J = 3, J = 4, and J = 6 rotational levels 
of H2 allowed us to plot the excitation diagram of H2 (see Fig. [2]). Assuming that all three 
levels are populated by thermal collisions (LTE), we estimated the excitation temperature 
of the observed gas. The upper limit on the column density of the J = 6 level strongly 
constrains the value of the column density of the J = 4 level. Indeed, the population of 
the J = 4 level should correspond to the excitation temperature given by the ratio of the 
column densities of the J = 3 and J = 6 levels. We assumed that the population of the 
J = 3 and J = 6 levels follow the Boltzmann law, which corresponds to a linear fit to the 
points on the excitation diagram (Fig. [21): 



where N{H2)j=i is the column density, gi is the statistical weight, and Ei is the energy 
of the J = i level; k is the Boltzmann constant and T^x is the temperature defined as the 



3.3. Physical properties of H2 




(2) 
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excitation temperature. Since we only have upper limits on the column densities of the 
J = 4 and J = 6 levels, only an upper limit on the excitation temperature is relevant. Since 
the temperature is inversely proportional to the slope on the excitation diagram, in order 
to obtain the upper limit on the excitation temperature, we considered the lower value of 
the column density of the J = 3 level (i.e., the measured value minus la) and the upper 
limit on the J = 6 population level (that procedure yields the minimum slope / maximum 
temperature and corresponds to the solid line on Fig. [2]). We thus find that the excitation 
temperature of H2 should be lower than 570 K. In this case, the column density of the J = 4 
level should be lower than 1.6x10^'^ cm~^ which is lower than the value obtained from the 
3a upper limit (see Table [2]). 

Under the assumption that the H2 emission is optically thin, that the emitting H2 is 
in LTE at a temperature of about 570 K, and that the source size is equal or smaller than 
VISIR^s beam size, we derived upper limits on the H2 mass (see Table |2]): 



M„„, = / X 1.76 X 10^2° 



•■gas 



{hc/A7vX)A^iXu{T) 



M, 



) 



(3) 



where F^i is the line flux, d the distance in pc to the star, a;„(T) is the fractional population 
of the level u at th e temperature T in LTE (for details on the calculation method, see 



van Dishoeck 



19921 ). / is the conversion factor required for deriving the total H2 gas 



mass from the H2-ortho or H2-para mass. Since Mj|/2=Mj|/2(ortho)+M//2(para), then / 
=l+l/(ortho/para) for the S(l) line (a H2-ortho transition) and / =l+ortho/para for 
the S(2) and S(4) lines (H2-para transitions). The equilib r ium o rtho-para ratio at the 



temperature T was computed using Eq.(l) from iTakahashil (120011 ). 
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4. Discussion 



Our previous high- resolution spectroscopic observations with VISIR of the S(l) pure 
rotational line of H2 at 17.035 fim of HD 97048 revealed the presence of significant amounts 
of warm gas in the surface layer of the disk in the inner 35 AU from the central star 
(Paper I). This detection confirmed that HD 97048 is a young object surrounded by a 
circumstellar disk at an early stage of evolution, still rich in gas. 

In the present paper, we report the non-detections of the S(2) and S(4) lines of H2 in 
the disk of the Herbig Ae star HD 97048. Using VISIR high-resolution spectra, we compute 
upper limits for the column densities of the J = 4 and J = 6 rotational levels. Using the 
results obtained from the analysis of the S(l) line (column density of the J = 3 level), and 
assuming that the J = 3 to J = 6 levels are excited by thermal collisions we derive an 
upper limit on the excitation temperature of H2 of about 570 K. 

This constraint on a temperature lower than 570 K for the gas allows us to set a 
stronger constraint on the mass of warm H2 (< 0.1 Mjup; 1 Mjup ~ 10~^ Mq) in the inner 
35 AU of the disk than that estimated in Paper I. Indeed, in Paper I, we derived masses of 
the warm gas in the range from 10~^ to nearly 1 Mjup, for temperatures in the range from 
150 K to 1000 K and assuming LT E. The mass o f warm gas we estimate is also lower than 



that we derived from the paper by 



Lagage et al. 



(120061 ) who estimated the mass of H2 in 



the disk to be of about 0.01 Mq. By assuming that t he surface density S follows a power 



law S(r) = S(r/370)'' with an index q equal to —3/2 (jLagage et al 



20061 ). we estimated in 



Paper I a minimum mass of gas in the inner 35 AU of the disk to be of the order of 3 Mjup. 
But it should be pointed out that mid-IR H2 lines are only probing warm gas located in the 
surface layer of the disk, when a higher mass of colder gas is expected to be present in the 
interior layers of the disk. 



However, from their disk model. 



Carmona et al 



(I2OO8I ) concluded that the expected 
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peak flux of tlie S(l) line at 17.035 /iin, observed at a spectral resolution of 20 000, should 
be less than 0.3% of that of the continuum at temperatures higher than 150 K, and thus 
should not be observable with the existing instrunaents. Those authors used a two-layer 



model (IChiang fc Goldreich 



1997 



DuUemond et al. 



20011 ) of a gas-rich disk (column density 



of N(H2)=10^'^ cm ^) seen face-on, located at 140 pc from the Sun, with LTE for the gas 
and dust, Tgas='^dust, and assumed a constant gas-to-dust mass ratio about 100. 

As shown in Paper I from the observation of the S(l) line in the disk of HD 97048, such 
a detection can only be e xplained if the physical conditions of the gas differ from those used 



by 



Carmona et al 



(120081 ). In Paper I, we proposed several ways to explain our detection. 
First, assuming equal dust and gas temperatures (LTE), we estimated gas-to-dust mass 
ratios much larger than the canonical value of 100. This hypothesis can not be ruled out 
by the non-detections of the S(2) and S(4) lines. One possible interpretation to explain a 
very high gas-to-dust ratio is that the dust is partially depleted from the uppermost disk 
surface layer, where the H2 emission originates. The spatial decoupling between the gas and 
the dust may be due to low densities in the surface layers, or dust settling and coagulation 
into larger particles. The physical conditions may thus rapidly differ from the LTE ones. 



i.e. Tgas > Tdust- In the upper disk surfac e layers, photoelectric heating can thus p 



significant role in the gas heating process (iKamp &: DuUemond 



2004 



Jonkheid et al. 



ay a 



2003). 



In addition, photoelectric heating will be much efficient on small grains, suc h as PAHs, 



which show significant features in the different observations of HD 97048 (e.g. 



2007 



Doucet et al. 



and reference therein). Our upper limit o n the excitation temp e ratur e is consistent 
with photoelectric heating. Indeed, as shown by lxamp fc DuUemond (2004), photoelectric 
heating is characterized by a temperature limit (~1000 K) since for high stellar UV flux, 
grain ionization can be so large that the photoelectric effect is less efficient and Tgas does 
not increase anymore. 
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However, our observations do not allow us to rule out the possibility of a gas heating 
due to other excitation mechanisms such as UV pumping or X-ray heating. Indeed, UV 
and X-ray heating are likely possibilities that can heat the gas to temperatures signicantly 



hotter than the dust. 



Nomura fc Millarl (120051 ) modeled UV heating of circumstellar disks 



around T Tauri stars and the resulting H2 emission. In their model the temperature can 
be of about 1000 K at 10 AU from the central star, and the predicted H2 S(2) line flux is 
about 2xl0~^^ ergss"^ cm~^. This value is lower than our upper limit on the S(2) line flux, 
but their model was for a less massive star and disk. 

The upper limit on the excitation temperature of about 570 K we foun d for H-;. in the 



disk o f HD 97048, could also be consistent with X-ray heating. Very recently. 



Ercolano et al. 



( I2OO8I ) demonstrated from their two-dimensional photoionization and dust radiative transfer 
models of a disk irradiated by X-rays from a T Tauri star, that the uppermost layers of 
gas in the disk could reach temperatures of 10^ K at small radii (<0.1 AU) and 10^ K at 
a distance of 1 AU. The gas temperatures decrease sharply with depth, but appear to be 
completely decoupled from dust temperat ures down to a column depth of ~5xl0^^ cm~^. 



These results are consistent with those of 



Glassgold et al. 



(I2OO7I ) who computed that at 



20 AU from the central star, the temperature can reach 3000 K at the surface of the disk 
before dropping to 500-2000 K in a transition zone and then to much cooler temperatures 
deep in the disk. By extrapolation, these models suggest that in the inner 35 AU from 
the central star, the upper layers of the disk could reach temperatures consistent with the 
upper limit of about 570 K we find for HD 97048. 

The kinetic temperature of the gas is generally supposed to be given by the population 
of the first rotational levels (namely J = and J = 1) since their critical densities are 
relatively low, and can not be higher than the excitation temperature given by the higher 
energy levels. Thus, in any case, even if the H2 we observe around HD 97048 is not 



-14- 



essentially excited by thermal collisions, the populations of the three levels derived from 
our observations give strong constraints on the gas kinetic temperature. Our upper limit 
on the temperature of about 570 K is thus reliable whatever the mechanisms responsible 
for the excitation of the observed gas. We stress that the excitation temperature we find 
for H2 from the VISIR observations of HP 97048 is clos e to that derived from the TEXES 



observations of AB Aur (670 K) by iBitner et al.l (120071 ). which reinforces the similarities 
between the two stars (see Paper I). 

In order to better constrain the physical conditions of H2, high-sensiti vity and high- 



resolu tion space-based spectrographs in the mid-IR would be required (e.g. 



Boulanger et al. 



20081 ). In particular, only the observation of the S(0) emission line near 28/im would allow 
us to clearly constrain both the kinetic temperature of the observed gas and the physical 
mechanisms responsible its excitation. 



This work is based on observations obtained at ESO/VLT (Paranal) with VISIR, 
programs' number 079.C-0839A and 079.C-0839B. CMZ is supported by a ONES fellowship 
This work was supported by Agence Nationale de la Recherche (ANR) of France through 
contract ANR-07-BLAN-0221. We also thank Programme National de Physique Stellaire 
(PNPS) of CNRS/INSU, France, for supporting part of this research. 
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Fig. 1.— VISIR spectra of HD 97048 at 17.035 //m {top panel), 12.278 //m {middle panel) 
and 8.025 //m {bottom panel). Left panel- continuum spectra of the asteroid and of the target 
before telluric correction. Central panel: full corrected spectra: dotted lines show spectral 
regions strongly affected by telluric features. Right panel: zoom of the region where the H2 
hues should be observed (dashed vertical hues) . The spectra were corrected neither for the 
radial velocity of the targets nor the Earth's rotation velocity. 
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Fig. 2. — Excitation diagram for H2 towards HD 97048. If the three rotational levels are populated 
by thermal collisions, their populations follow the Boltzmann law, and the upper limit on the excitation 
temperature is given by considering the lower limit on the column density of the J = 3 level and the upper 
limit on the J = 6 population level (solid line). Thus, in this case, the gas temperature should be lower than 
570 K. 
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Table 1: Summary of the observations. The airmass and seeing intervals are given from the 

beginning to the end of the oljsorvations. 
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Table 2: Fluxes, intensities, luminosities and column densities of each observed mid-IR tran- 
sitions of H2. The mass of H2 is calculated for a temperature of about 570 K. Xobs is the 
wavelength of the centroid of the line. Vup and Jup are respectively the vibrational and 
rotational upper levels of the transition of interest. 
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■^i The error bars on the continuum fluxes are 3a error bars. 



